Generalized Parton Distributions and Color Transparency 

Phenomena 

S. Liuti 1 '! and S. K. Taneja 1 '! 

1 University of Virginia, Charlottesville, Virginia 22904-4114, USA. 

Abstract 

We study the structure of generalized parton distributions in impact parameter space with the 
aim of determining the size and role of small transverse separations components in the quarks 
wave function. We analyze the relation between transverse momentum components and transverse 
separations. Wave functions with large transverse momentum components can simultaneously 
reproduce the behavior of the Dirac form factor at large momentum transfer, and of the deep 
inelastic structure functions at Bjorken x —* 1. The presence of large momentum components 
does not ensure, however, the dominance of small transverse distances at large x. We suggest that 
experiments measuring the attenuation of hadrons in the nuclear medium, or the onset of color 
transparency, can provide an alternative source of information on generalized parton distributions, 
by mapping out the behavior of the transverse components of the wave function. 
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I. INTRODUCTION 



A most intensively studied question in Quantum ChromoDynamics (QCD) is the space- 
time structure of high energy exclusive reactions. In perturbative QCD (pQCD) approaches 
the hypothesis is made that elastic processes are dominated by the Fock space components 
of the hadronic wave function with the minimum number of quarks (anti-quarks). In the 
hard scattering picture these are assumed to be located within a small relative transverse 
distance < 1/ y/Q 2 , Q 2 being the (high) four-momentum transfer squared in the reaction 0. 
After the hard scattering takes place, a normal hadron is recovered. The validation through 
experiments of this picture - the dominance of the short separation components - is however 
not straightforward, as witnessed for instance bv the recent controversial interpretations of 
experiments sensitive to helicity selection rules 0, |j| . 

It has been long surmised jj], 0] that performing (quasi)-elastic reactions off nuclear 
targets might ease this impasse. Nuclei can in fact function as both "passive" or "active" 
probes for the small partonic separation components. Small distances can be filtered at 
finite Q 2 , by impeding the passage of large separations because of the strong interactions 
occurring in a nucleus with A >> 1. In this case nuclei play an active role (if. One can also 
simply study the passage through the nuclear medium of the small separation component 
selected at asymptotically high Q 2 , namely the phenomenon of Color Transparency (CT). 
CT is attained in principle when the ratio, Ta(Q 2 ) = oai 1 Aoh-, of the nuclear cross section 
to the one for scattering in free space, reaches the value of 1 at some given high value 
of Q 2 . From a practical point of view, current searches for CT might appear to be in a 
stall as all experiments performed so far seem not to show either any systematics or any 
marked trend for the onset of this phenomenon 0]. However, at the same time, none of 
the experiments does contradict current predictions for CT. A large body of theoretical 
work exists, that generated an active field in the last decade (see 0, || 0] for extensive 
reviews). New observables have been proposed since the original measurements, such as the 
attenuation cross section [Hf , the systematics of oscillations in pion production reactions , 



or the sensitivity to the correlation length in exclusive vector meson production . Based 



on these new hypotheses, more specific searches have been conducted, which in some cases 
seem to better identify the phenomenon [l2l . Il3l FLU Critical observations about the 
key assumptions underlying the onset of CT have also recently emerged, a most disquieting 
one being the suggestion that the transverse size of the exclusive process might not even be 
small due to a "suppression of the Sudakov suppression" jT(| . 

Whether or not a pQCD description of the hadrons holds at the Q 2 values presently 
available, it has now become imperative to investigate the basic question of the existence 
and observability of small size hadronic configurations. In this paper we use as a theoretical 
tool Generalized Parton Distributions (GPDs) that, as pointed out recently [T3|, can give 
us direct access to the transverse spatial extension of the hadronic components. GPDs can, 
in fact, be shown to be diagonal in the impact parameter representation, and they can 
therefore be interpreted as - Fourier transformed in the transverse plane - joint probability 
distributions of finding a parton with longitudinal momentum fraction x, and transverse 



separation from the proton's center, b [17|. The original result, obtained for £ = 0, £ being 
the longitudinal momentum transfer fraction, was shown to hold also in the £ 7^ case in 
[l8l | (see also 3| for a laboratory frame based description). 

The new insight provided by the impact parameter space interpretation allows one to 
make a connection between GPDs and the phenomena of CT and nuclear filtering. Within 
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this new framework, one can in fact address directly the question of what transverse distances 
dominate the hard exclusive scattering cross section. So far, a quantitative knowledge of the 
size of the hadronic configurations involved in a reaction has been in fact very elusive, or at 
the best, loosely inferred from experiments. 

This paper is organized as follows: in Section |H] we define GPDs and show their connec- 
tion, at leading order in Q 2 , to both the impact parameter dependent parton distributions, 
and to the intrinsic transverse momentum distributions. In Section imi we evaluate the role of 
large transverse momentum components and, conversely, of small interparton separations. 
A few quantities are calculated, namely the GPDs at £ = 0, the Dirac form factor, and 
the the deep inelastic structure function. For each of these quantities the impact of differ- 
ent behaviors in the partons wave functions in both transverse momentum, and transverse 
distance, is evaluated. In Section HVl we calculate the ratio Ta{Q 2 ) for electron nucleus scat- 
tering experiments in terms of GPDs, and we suggest a method for extracting the transverse 
interparton distance distribution. Finally, we draw our conclusions in Section 



II. SOME OLD AND NEW DEFINITIONS 



GPDs were initially introduced |20l . |21| in order to better describe in a partonic language 
the concept of orbital angular momentum carried by the nucleon's constituents. They can be 
accessed experimentally in exclusive hard scattering experiments. The most straightforward 
way has been to date Deeply Virtual Compton Scattering (DVCS), that has become the 
prototype for illustrating different aspects of GPDs. There exist a few reviews to date on 
this subject 0, H3- I n this Section we summarize the properties of GPDs that are relevant 
for our discussion of exclusive scattering from nuclear targets, as e.g. probed in the reaction: 
eA — > e'p(A — 1). In particular, we define: 

• The Impact Parameter dependent Parton Distribution Function (IPPDF) and its re- 
lation to the observables, F\, the Dirac form factor, 1/W2, the DIS structure function, 
and H, the "non-spin flip" GPD. 

• The rms impact parameter, interparton separation, and radius for the proton, in terms 
of both IPPDFs and GPDs. 

• The intrinsic transverse momentum distribution, and its relation to both the IPPDFs, 
and to Fx, 1/W2, and H. 

The latter issue opens a new perspective on GPDs as "tools" to study transversity. 



A. Generalities 

We consider the leading order in Q 2 (twist-2), where one can define four universal GPDs, 
namely H, E and H, E. The latter two are related to the quarks helicity distributions. In 
this paper we focus on the unpolarized case described by H, E. The relevant kinematical 
variables are: P and P', the initial and final nucleon's momenta in the exclusive process, 
P = (P + P')/2, the average nucleon momentum, k the active quark momentum, and q the 
virtual photon momentum (Fig.l). While the (upper) hard scattering part of the diagram 
varies in principle for different reactions, the information necessary to describe GPDs is 
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contained in the flower) soft part (the different types of matrix elements related to GPDs 
are described e.g. in [22]). 1 One can define four independent invariants using these 
variables: 2 x = k + /P + , £ = -A+/2P+, Q 2 = -q 2 , and t = A 2 , where A = P - P', is the 
four-momentum transfer in the process. Differently from inclusive Deep Inelastic Scattering 
(DIS), the four-momentum transfer, A, does not generally coincide with the virtual photon 
four-momentum, q. By properly devising experimental measurements that reproduce this 
situation, namely q ^ A, one has therefore the unique possibility of "zooming" in on partonic 
components, due to the high resolution provided by large Q 2 (in the multi-GeV region), at 
the same time scanning the hadronic target's spatial distribution, through the dependence 
on the variable A. We denote by Hi and Ei the GPDs for a quark of type i. Hi and Ei 
satisfy the following relations: 

j dxHi(x,£,t) = Fi(t), (la) 
J dxEi(x,£,t) = F^t), (lb) 



independently from £ (see j22j,|23|). F^ are the Dirac and Pauli form factors for the quarks 
of type i in the nucleon. They are readily obtained from the measured proton and neutron 
form factors by using isospin symmetry |22j. For x > 0, 

fl-<(x,0,0)=ft(x), (2) 

i.e. one recovers the Parton Distribution Functions (PDFs), qi(x), in the in the forward 
limit. E does not have an equivalent connection to parton distributions. 



B. Impact Parameter Space 



GPDs can be related jl7( to the Impact Parameter space dependent PDFs (IPPDFs) (see 
e.g. Ref . [iif ) . An IPPDF is defined as the joint distribution: drii/dxdh = qi(x, b) - the 
number of partons of type q with momentum fraction x = k + /P + , located at a transverse 
distance b (b is the impact parameter) from the center of P + of the system. 3 The latter is 

defined as R = (S^ ^t^n) / P + , where x% = (x°; x = b, x 3 ) is the space-time variable for the 

n — th parton, and k£ = P + ■ The usage of Light Cone (LC) variables is motivated by 

the fact that by viewing x + as the time variable, P + ,p+ as the total and partonic masses, 
respectively, and R as the center of mass of the system, P + R is the generator of longitudinal 
boosts that leave the physics on the x + = surface invariant. An intuitive connection 
with non-relativistic logics is naturally provided. The interpretation of observables, e.g. 
structure functions, in terms of parton density distributions relies on the validity of the 
Impulse Approximation (IA) whereby the hard scattering happens instantaneously in x + . 
We work in this hypothesis by assuming that it is satisfied in the kinematical limit for the 



1 We do not take into consideration initial state nuclear medium modifications. These are part of a more 
general problem, and they will be addressed in a forthcoming paper |24j 

2 We use the notation: a M = (a ; aj_, a 3 ), and a ± = (a ± a 3 ) 

3 All transverse variables are denoted here by bold-face characters. 
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coherence length of the virtual photon given by: i + « i k XjlM^xsy This amounts 
to considering only large values of xbj- At low xbj, one has to consider both IA-type 
scattering and of processes including photon splitting into quark-antiquark pairs before the 
nucleon where it has been recently shown that FSI might hamper the interpretation of 
the deep inelastic measurements in terms of probability distributions [2^]. This is part 
of a wider set of problems that transcend the goals of this paper. In particular, it will 
be important to study in the future whether parton interactions (or off-shellness) affect the 
transverse distributions that we address in this paper, as this might complicate the procedure 
of extracting holographic images of the proton as proposed in jlil l27l|. 

By definition |25[, The DIS structure function, uW2(x), and the form factor, Fi(Q 2 ), can 
be written in terms of the IPPDFs respectively, as: 

uW 2 {x) = ^e 2 x I d 2 b qi (x,b) (3a) 

i J 

Fi(A) = J2 ei J rf2beiA b / dx qi (x,h), (3b) 

where e; is the quark's charge, and we have disregarded the Q 2 dependence of vW 2 - Notice 
that A is the conjugate variable of b. 

1. b distribution 

In ^?|]; the observation was made that for a purely transverse four momentum transfer, 
namely for A = (A = 0; A, A 3 = 0) and £ = 0, H^x, 0, —A 2 ), and qi(x, b) can be related 
as follows: 

tf ( (,,0,-A') = /A^(,,b) (4b) 

Since qi(x, b) satisfies positivity constraints and it can be interpreted as a probability dis- 
tribution, Hi(x, 0, —A 2 ) is also interpreted as a probability distribution, namely the Fourier 
transformed joint probability distribution of finding a parton i in the proton with longi- 
tudinal momentum fraction x, at the transverse position b, with respect to the center of 
momentum of the nucleon. Since in what follows we will consider only the case £ = 0, we 
define: H^x, 0, -A 2 ) = H^x, A). 

The root-mean-square (rms) distance of the configuration with momentum fraction x, 
from the center of the proton is defined as the square root of: 

( b 2 (x) ) = N J d 2 hq(x,h)h 2 , (5) 
M being a normalization factor, i.e. Af~ l = J d 2 h q(x, b). Similar relations to the ones in 



Eqs. fl4al4b|) . were obtained for £ 7^ in [18|. An alternative approach in terms of Wigner 
distributions, formulated in a covariant model that allows one in principle to perform calcu- 
lations in the laboratory frame was given in The role of the variable £ in the context 
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of a covariant approach, is an important subject per se, especially in reference to Initial 
State Interactions (ISI) in nuclei, since it impacts the particles' off-shellness. We address 
this subject in detail in 



2. y distribution and hadronic radius 

One can also define the Fourier transform with respect to the variable y = b/(l — x), 
the transverse position of parton i with respect to the center of momentum of the system of 
spectator partons |25| : 

«,(,,£) = j /W^,b) (6) 

7i can be obtained from H by rescaling the functions arguments as follows: 

H i (x,A)=H i (x,A(l-x)) (7) 

Since y measures the interparton separation, it is useful to define also the inverse Fourier 
transform: 

H t (x,y) = j ' dk!e-^Hi{x,V) = q{x,y(l-x)). (8) 

The variables k = A(l — x) and y are Fourier conjugates, as well as A and b = y(l — x). 
7i(x, y) is the probability of finding a parton with momentum fraction x, at a transverse 
distance y from the center of momentum of the spectator partons. The rms transverse 
distance between partonic configurations with momentum fractions x is obtained as the 
square root of: 

(y 2 (x))=M J d 2 bg(x,b(y))y 2 , (9) 

i.e. from the same observable H used to obtain (b 2 ). Af is defined analogously to Eq.(|5|). 
The radius of the system of partons, which is the useful quantity for our CT study, is: 

(r 2 (x) ) 1 / 2 = MAX {(b 2 ^)) 1 / 2 , (b 2 ^)) 1 / 2 ^} (10) 



This definition coincides with the one given in [l2j for x — > 1. In FigElwe sketch a cartoon of 
the different transverse distances discussed so far. The figure shows three different scenarios 
for the behavior of such distances with increasing x, corresponding to the models discussed 
in the following Sections. 



3. Intrinsic transverse momentum distribution and spectral function 

The function 7ii(x, k) defined in the previous subsection, although being related through 
a Fourier transform to the interparton separation, y, is conceptually different from the 
"intrinsic transverse momentum" distribution. 

Although we started from writing the structure function and the form factor in terms of 
IPPDFs (Eqs.([3a|) and (|3b|). respectively), their definition is usually given in terms of a LC 
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wave function, k) as [281 ] : 

vW 2 [x) = J2 e2 i x I d2k f( x ' k ) ( lla ) 

i J 

F(A) = Yl ei J d * k f Q <*z/(x,k,k+(l -x)A), (lib) 

where: 

/(x,k,k') = 0*(x,k)0(x,k'), (12) 

with k' = k + (1 — x)A, is a non-diagonal intrinsic momentum distribution. The diagonal 
term can be written as: 

/(x,k)=|0(x,k)| 2 . (13) 

By comparing Eqs. ljllal II lb)) , with Eqs. ()3al I3b|) . we find that the relation between the 
transverse momentum and the transverse separation of quarks inside a hadron is obtained 
through a non-diagonal distribution in transverse coordinate space, q(x,b,h'): 

f(x, k) = J d 2 b J d 2 b' e ik - (b - b ' } q(x, b, b'), (14) 

where we define: 

q(x,b,b') = **(z,b')*(z,b), (15a) 

q(x,b,b) = | i&(x, b) | 2 = q(x, b), (15b) 

\l/(x, b) being a wave function in the mixed x and b representation, related to the LC wave 
function by a transverse Fourier transform. This result can be also implemented in studies 
of processes directly sensitive to k. It is useful for this purpose to consider the average 
transverse momentum squared, defined as: 

( k 2 (x) ) = N J d 2 k f(x, k) k 2 , (16) 

with AT- = /d 2 k /(.. k). In Section Owe present an initial exploratory study of the 

relative roles of the intrinsic transverse momentum and the transverse separation using 
different types of wave functions. 

Finally, it should be mentioned that Eq. ljllaD corresponds to a projection obtained from 
the partons covariant spectral function, S(k) |29l 30l 31]: 



vW 2 {x) = Y,e 2 x J d 4 k S{k)s(x-^j 

= J2 e " x j d ' k J dP xS(k), (17) 

where the integration over k~ in the last line is replaced by an integration over the invariant 
mass of the spectator system, P\ = (P — k) 2 (k~ and P\ are related by). The main reason 
for introducing the spectral function along with the intrinsic momentum distributions is 
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that it encodes information about the energy and momentum of the spectator system. The 
spectral function can be defined in fact as the joint probability distribution of finding a 
parton with LC components k + and k, leaving the spectator system with a mass M\ = P\. 
The usage of the spectral function is most important in problems where the particles' off- 
shellness/transversity plays a role, a typical example being deep inelastic scattering from 
nuclei |24l |29l I30L l31|. In the following Sections we show a model of GPDs that implements 
a function S(k) within relativistic IA. 



III. MODELS FOR IMPACT PARAMETER DEPENDENT PARTON DISTRIBU- 
TION FUNCTIONS 

We evaluate the functions defined in the previous Section at leading order in Q 2 , using a 
two component, or spectator, model (Fig. lb). The main assumptions of the model are: i) the 
validity of the hand-bag diagram; ii) a description of dynamics in terms of a nucleon- quark- 
spectator system vertex. The minimal Fock component is described by a quark-diquark 
vertex function. We describe the active quark as off-shell and the diquark as on-shell, with 
mass Mx = yPx- 

The spectator model is both widely used and well tested in hadronic physics. In particular, 
in (2^1 it was used to describe DIS in nuclei; in |3(| parametrizations for both the polarized 
and unpolarized proton and neutron structure functions were given; more recently it was 
applied to DVCS [H, H^j. Our usage of the model other than for practical reasons, is 
motivated by the clear handling of the interplay between the particles' off-shellness and 
their energy and transverse momentum. This turns out to be an important property when 
describing the deep inelastic structure of nuclei j24|. By considering the components in 
Fig. lb, one has 



k 



.12 



k + k~ - 
(xP + ) 

k' + k'~ 
(xP + ) 



k 2 = 

M 2 
-fk- 



M x 



(1 - x)P* 

A) 2 = 



M 2 + A 2 M 2 + k 2 



P+ 



(l-x)P+ 



(k+A) 



;i8a) 



;isb) 



where A 2 = —t; M and M x are the proton and the diquark masses, respectively. Hence, 
for the denominators in Fig. lb: 



D(x,k) 



M x x - 



x 



(19) 



where we defined M. 2 X = M 2 x — M x /(1 — x) — m 2 /x, m being the struck quark mass. We 
consider the remnant diquark system to be a scalar. The form of the vertex function depends 
on the assumptions made in treating the quark's spin. Since we are ultimately interested in 
determining the radius of the proton as extracted by using the spin averaged quantity H, 
we assume that the vertex function is a scalar in Dirac space described by a function g(k 2 ). 
The two-component wave function becomes: 



(x,k) 



g{k 2 ) 

D(x,k) 



(20) 
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Notice that, because of the relation between the struck parton's virtuality and transverse 
momentum (Eq. ([18bp ). one can write as a function of the invariants: k 2 , x, M x . One can 
therefore identify the relativistic spectral function in Eq. (J17|) with: 

S(k)=\<f>(k 2 ,M x ,x) | 2 . (21) 



A. GPD, Structure Function, and Form Factor 

The GPD, the structure function, and the form factor can be evaluated in terms of the 
function in Eq. (j2U|) . according to the definitions in Section Hi B 31 One has: 

H(x,A) = y"c/ 2 k0*(x,k)0(x,k+ (1 -x)A) (22a) 
uW 2 (x) = J d 2 k | 4>(x,k) | 2 (22b) 



Fi(A 2 ) = j d 2 k J dx4>*(x,k)4>(x,k+ (1 -x)A), (22c) 

Notice that we incorporated in the function <f>, a factor \/x/(l — x), that accounts for both 
the initial flux, and the phase space for the process in Fig.l. Moreover, we drop the index i 
in H, since we fit our functions directly to the proton valence structure function and form 
factor, without distinguishing between u and d quarks contributions. For the vertex function 
we used the form: 

where g is fixed by the normalization, and A is a cut-off parameter. The values of A, m, 
and Mx are determined by fitting both the large four-momentum behavior of the form 
factor, and the shape of the DIS structure function at large Q 2 . The latter is determined 
in our model calculation at very low Q 2 (Q 2 < 1 GeV 2 ), where the structure function is 
given only by its valence component, the sea quarks and gluons being generated through 
perturbative evolution (see e.g. j^)- The values of the parameters are: m = 0.3 GeV, 
Mx = 1.100 GeV, A = 0.73 GeV. In Figures 3, 4, 5, and 6 we present our results obtained 
using Eas. (l^0l2^al22bl2^cli3l) . 

Fig.3 shows the quantity: | $(k) | 2 , where 

$(k) = f dx(f>(x,k). (24) 
Jo 

Also shown for comparison are results obtained using the parametrizations from Refs.|3 

t l37T |. The wavefunction in Eq. (J20|) . has a similar behavior to the parametrization used in 
, in that it displays high momentum components with a oc 1/k 4 behavior, although in 
a slightly larger amount. In order to quantify the "hardness" of the distributions, in Table 
|l]we show the average transverse momentum values obtained in each model. In the Table 
we also display the contributions to (k 2 ), from the large momentum components calculated 
by setting the lower limit of integration in | k | equal to 1 GeV and to 2 GeV, respectively. 
These values reflect the behavior of the curves in FigOJl in that the distribution in Ref. |21 



is truly of the "soft" type; the ones in Ref. [36| and in this paper are characterized by a long 
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TABLE I: Contribution of large transverse momentum components from different GPD models: 
values of (k 2 ) (first column); percentage contribution to (k 2 ) from k > 1 GeV components (second 
column), and from k > 2 GeV components (third column). 





((k 2 ) 1 / 2 MeV) 


% value for k min = 1 (GeV) 


% value for k min = 2 (GeV) 


This paper 


360 


14 


1.5 


Radyushkin '98 [35] 


290 


0.1 





Stoler '02 [36] 


286 


2.0 


1.0 


Burkardt '04 [22] 


415 


15 


6 x 10~ 2 



high momentum "tail" ; finally, we show results obtained with parametrizations of the type 
proposed in Ref.j3?]], where by accounting for gluonic interactions through the introduction 
of an effective gluonic mass in the initial two component model of (35[, a further suppression 
in the large x behavior is found, of the type ~ exp [(1 - x) n A 2 ] i with n > 2. We find that for 
n = 2, 3 the distributions obtained following the idea of Ref.|37[, are neither entirely "soft" 
nor "hard" , but they display larger "intermediate" momentum values, hence they give rise 
to a large (k) 2 . 

The function H (x, A) is shown in Fig. 4, plotted vs. A for different values of x. A small 
(x = 0.07), an intermediate (x = 0.36) and a large (x = 0.88) value of x are indicated in 
the figure, from which one can also deduce that the value of H(x, 0) traces the expected 
behavior for the valence structure function (namely the curve at x = 0.07 has a much lower 
value for A —* than the curve at x = 0.36). Values of x > 0.9 decrease one more order of 
magnitude and they are not shown in the figure. 

In Fig. 5 we show the form factor .F\(A 2 ). This is compared to calculations using the 
models of Ref.|2l[, Ref.Q, and to the dipole model: F^A 2 ) = + A 2 /0.71) 2 , with 

T = (A 2 /4M 2 /x P + l)/(l + A 2 /4M 2 ). We find, in agreement also with Ref.@, that the large 
momentum components of the wave function are responsible for the large A 2 behavior of 
the form factor. Such components are effectively comparable to the hard contribution from 
pQCD calculations jjj^]. Although they constitute a small fraction of the average transverse 
momentum as shown in Table I, the contribution of large k is much larger in the nucleon form 
factor. This is shown in Fig.6, where we plotted the ratio: i*i(A, k max ) / Fi(A, k max = oo). 
Only w 70% of the form factor is reproduced by stopping the integration at k max = 0.7 GeV, 
by using a wave function with hard components, whereas the ratio shown in the figure is 
already saturated in models based on a Gaussian behavior. 

Similar conclusions can be drawn by investigating the large x behavior of the deep inelastic 



structure function, Eq. (j22b 



B. Impact Parameter Space Parton Distributions 

Next, we study how the behavior of the wave function with transverse momentum, k, 
reflects on the transverse space distribution, q(x, b). This is obtained as a Fourier transform 
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of the function H(x, A) in the transverse plane: 4 



1 r°° 

q(x, b) = — dAA J (bA) H(x, A) (25) 
2tt J 

where A = | A | and b =| b |, and Jo is a zero-th order Bessel function of the first kind. In 
order to evaluate q(x, b) with adequate accuracy, we parametrized H (x, A) as: 

A ) = t4^ TTT^- (26) 

(p 2 (x)+p 3 (x)A 2 +p 4 (x)A) 2 v ; 

The coefficients pi, i — 1,4 are given in the Appendix. The Hankel transform of /J has the 
following analytic form: 

q(x, b) = \ b Kl (^p 2 (x)/p 3 (x)b), (27) 

where K\ is a modified Bessel function of the second kind, and we set p 4 = (a more 
complicated, though similar, analytical form is obtained for p 4 ^ 0, which we do not display 
for ease of presentation). q(x, b) is shown in Fig. 7, plotted as a function of b for varying x. 

Using Eq.([27|) one can calculate the radius, the interquark separation, and the rms dis- 
tance from the center of the hadron, Eqs. (j5l9ll0|) . These quantities are shown in Fig.8. In 
Fig. 8(a) results for (fc 2 ) 1//2 from several models for the GPDs are compared. The difference in 
slope as x — > 1 is clearly visible. In Fig.8(b) we show our results for (b 2 ) 1 ^ 2 , (y 2 ) 1 ^ 2 , (r 2 ) 1 / 2 . 
One can see that the both rms radius, and the interparton separation do not vanish in the 
x — > 1 limit; rather, they tend to a finite (small) value of « 2.5 GeV" 1 (~ 0.5 fm). In Fig. 9 
we plot the average transverse momentum squared of a quark in the proton, Eq.([16|). There 
is a marked difference in the behavior of this quantity for the "soft" and "hard" type wave 
functions. The wave function proposed in j3?J gives manifestly divergent results for large x, 
although, as one can see from Table I, the integrated value of (k 2 ) is finite. Our results are 
consistent with previous findings using the diquark/spectator model 30l . 



What components of the wave function are responsible for this behavior? In Fig. 10 we 
give a partial answer to this question by plotting the contribution to both b (upper panel) 
and r (lower panel) of the short distance components. This is obtained by stopping the 
integration over b in the definition of the rms quantities at decreasing values of the upper 
limit: 2, 1, and 0.5 GeV -1 , respectively. It is clear that the short distance part of the 
wave function (or the large k components) largely determine the size of the radius at large 
x. This is particularly evident from the lower panel of Fig. 10, where results are presented 
using the ratio: (r 2 (b max )) / (r 2 (oo)) . While models including a hard component saturate 
this ratio for b < 0.5 GeV -1 already at x — 0.8, models using a "soft" wave function of the 
type of Ref. [35| do not show a tendency to saturation until very large x. These results open 
up a somewhat puzzling situation: In order to describe the form factor at large Q 2 , one 
needs high k components in the wave function. This is both a prediction of pQCD based 
calculations 0,01; an d it can be also shown phenomenologically (Fig. 5) (see also Ref.|36j]). 
High k components are also responsible for a decrease in the hadronic transverse size as x 
increases (Figs. 8, 9 and 10). However, one can envisage a model such as the one utilized in 



For a circularly symmetric function of A this is a Hankel transform 
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this paper that has the amount of high momentum components necessary to reproduce the 
Dirac form factor, and that at the same time, predicts a small but not vanishing transverse 
size for hadronic configurations at large x. 

In concluding this Section, we reiterate that the usage of GPDs for determining the 
spatial structure of hadrons through q(x,b), has been so far largely theoretical, since these 
quantities are only indirectly deduced from a variety of delicate measurements ([H,|3^]). In 
the following we suggest that an additional source of information could be given by studies 
of Color Transparency and related phenomena such as Nuclear Filtering. These are in fact 
directly sensitive to the size of hadronic components. 



IV. EXCLUSIVE SCATTERING IN THE NUCLEAR MEDIUM 

In the previous Sections we discussed the relative importance of the soft and hard com- 
ponents of the light cone wave functions up to four momentum transfer Q 2 = 25 — 30 GeV 2 . 
Plausible scenarios confirm the idea that ~ 30% of the form factor is given by the hard com- 
ponent of the wave function at four momentum transfer Q 2 = 25 — 30 GeV 2 (for consistency 
with the general notation we switch from A 2 to Q 2 in this Section). This result is indeed 
correlated with the dominance of small transverse sizes of hadronic configurations through a 
Fourier transform, however not straightforwardly as demonstrated by the different behaviors 
of the GPDs with x and A. As shown in Fig.ll where the average value of x, given by: 

* * \ L dxxH(x, A) 

aw A = ^ )lJ — L , 28 

1 ; J dxH(x, A) 1 ; 



is plotted vs. A = yQ^, the model in Ref.|3J| is governed at large Q 2 by small b as well as 
by large x components (although its prediction for the form factor is substantially lower than 
the data, as shown in FigEl and the radius predicted for hadronic configurations diverges at 
large x). The approach based on the diquark model yields a form factor that is dominated 
by small b components but not exclusively by x — > 1 configurations. This behavior which is 
further clarified by the comparison of H (x, A) for the two cases in Fig. 12, is at the origin of 
the rather flat dependence with x of the transverse radius shown in Fig. 8. Also, it appears 
that the model of Ref . [37^ . approximately including the effect of a gluonic interaction, fixes 
the problem of an unphysically large radius while keeping a large x behavior as the one in 
Ref.|25j, at the expense of introducing a large k in the wavefunction (FigEJ. 

In summary, while it can be challenging to unambiguously disentangle the amount and 
nature of hard components responsible for the large Q 2 behavior of the hadronic form fac- 
tors, one might gain a better insight by requiring models to simultaneously describe the 
hadrons transverse spatial distribution, and in particular the possible onset of small trans- 
verse configurations. The diquark model presented in this paper seems to provide a better 
description than other models similarly based on the handbag diagram. 

If configurations with small radii indeed exist, they can be isolated in principle by per- 
forming CT and/or nuclear filtering type experiments. CT has been so far investigated 
in a variet y o f reactions including proton-hadron, 7-hadron, and electron-proton scattering 
T3I . 14j . Here we consider exclusive electron-proton scattering, where oscillations due 



to the interference of perturbative and non-perturbative contributions possibly present in 
the proton-proton and 7-proton exclusive cross sections, do not occur. Our result is however 
quite general and can be easily extended to other type of reactions. 
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Using the new perspective of GPDs in impact parameter space, we can now take a 
radical turn and devise a "minimal" approach that reduces the (nuclear) model dependence 
of the problem (an approach complementary to ours was recently suggested in 4y| , where 



the connection between CT and GPDs was explored in relation to diffractive dissociation 
of pions). As widely underlined in the literature, the prerequisites for the onset of CT, 
assuming that the scattering takes place off a small size configuration, are that: i) the 
interaction of the hadron in the nuclear medium is reduced due to the decreased gluon 
radiation from a small transverse size color dipole ji^, 41 1; ii) the small size configuration, 



not being stationary, will in principle evolve with time into a larger configuration. It however 
remains sufficiently small during the time it crosses the nucleus. Point ii) has probably been 
the most controversial one, although recent criticism to i) has also emerged ^6[. Another 
point, not extensively addressed in the literature (see however 0), is that the hadron suffers 
Initial State Interactions (ISI), or that in other words, the form factor for an off-shell nucleon 
is in principle different from the on-shell one. 

We can now break down the problem in several steps: We first assume that the scattering 
happens with an "unmodified" proton in the nuclear medium. We introduce a nuclear filter 
for the large transverse size components, by multiplying the IPPDF, q(x,b), by a square 
function: 

n m = / 1 b < b max{A) 

K} I b > b max {A) 
where b max (A) defines the size of the filter. This affects the GPD as: 

rb m ax(A) 

H A (Q 2 ) = / dbbq(x,b)J (bA), (29) 

Jo 

with A = s/Q 2 . The transparency ratio for (e, e'p) type reactions becomes 

" i 1 2 

L dxH A (x } A) 

T a(Q 2 ) = T2-' (3°) 

/ dxH(x, A) 

where H defined as in Eq. (j4b|) represents scattering in free space, i.e. it is calculated with 
b-max — oo. Based on this result, one can fit the available data, using different distributions 
q(x, b), and varying the parameter b max . In FigiT3]we show a few possible scenarios based on 
a simplified analytical model for q(x, b). We take q(x, b) = A(x) expf— a(x) b), where a(x) is 
taken with two different dependencies: (a) a oc (1 — x)/x, as in Ref . |35j] : (b) a oc (1 — x) 2 , 
as in Ref.jl3|. The effect of the filter is to produce both damping and oscillations in Ha- 
In FigfTBI we show for illustration, the ratio R = Ha(x, A)/H(x, A) plotted vs. A for two 
different values of x, in case (a) and (b), and for different values of the filter size. Our 
analytical calculation displays explicitly the damping and oscillations as: 

R = 1 - exp(-a b max ) (a b max J (Ab max ) + cos(Ab max )) . (31) 

Although this result does not correspond to a completely realistic situation, it allows us to 
understand for varying x, the different effects due to variations in the size of a, which in 
turn is a feature of GPDs from different models. Oscillations are clearly more pronounced 
at low x. When b max is large, of the order of 5 GeV -1 (~ 1 fm), transparency is attained. 
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When b max is small (1 GeV" 1 , ~ 0.2 fm), the medium is no longer transparent to the type of 
distributions considered, although distribution (a) giving a larger radius, is more suppressed 
than (b) at low A. 

The results of the first step of the analysis would give joint information on possible sizes 
of the filter, defined as y max = b max / (1 — x), and on the behavior of the IPPDFs, i.e. on the 
radius of hadronic configurations, (r 2 ) 1 / 2 . 

This procedure, systematically applied to a sufficiently large body of data including both 
existing and planned measurements 0, 0, Q, 4^], gives a much more direct test of the 



transverse sizes involved. Once this basic information is known from measurements - and 
not inferred from theoretical scenarios - one would be able to introduce more sophisticated 
modeling of e.g. rescattering and ISI. A somewhat similar point of view was taken in a 
pre-GPD context in Ref.|43| and applied to the data in 0,0 , where fits were performed 
based on a parameter, p = nARA&eff, with via the nuclear density, Ra the nuclear radius, 
and a e ff the effective hadron-hadron interaction cross section. 



V. CONCLUSIONS AND OUTLOOK 

We presented a study using the new theoretical insight provided by GPDs, of the in- 
terplay between the transverse variables in both momentum and coordinate space, k and 
b respectively, and of the longitudinal momentum fraction, x, in the proton form factor, 
Fx. Our study was aimed at establishing what type of components - hard vs. soft, small 
transverse distance vs. large - dominate the form factor based on different hypotheses for 
the GPDs used to describe Fx- 

We examined three types of distributions (Figs|3]and|H): 



i) Parametrizations based on Gaussian wave functions with argument oc (1 — x)A 2 |3j 

ii) Parametrizations based on Gaussian wave functions with argument oc (1 — x) n A 2 , with 
n > 2, 0; 

Hi) GPDs obtained using a diquark/spectator model consistent with a 1/k 4 asymptotic 
dependence of the wavefunction. 

We confirm the result initially obtained in Ref . |36| (Fig|SJ) that Fx at A > 2 GeV can be 
reproduced only by allowing for a sufficiently large amount of large k (hard) components. 
Cases ii), and Hi) satisfy this condition. Furthermore, we find that, for values of A 2 pa 25—30 
GeV 2 , the behavior predicted for the form factor fFigsIUl I§1 ITTj) is governed by: 



i) Large x (x ave (A = 5 GeV) « 0.75), k < 0.7 GeV in |3J 

ii) x — > 0, large k in jj^l; 
Hi) x pa 0.25, k < 2 GeV in our diquark/spectator model. 

The rms value of the proton radius, given by (r 2 ) 1/2 in Eq.jHIJ), was calculated by a 
two-dimensional Fourier transform of the GPD, i.e. using the function q(x,b) (FigJTJ). We 
pointed out that (r) 1 / 2 — > (y) 1 ^ 2 for x — > 1, where y is the variable defining the radius 
in Ref. |37| . and (r) 1 / 2 — > (b) 1 / 2 ^ (y) 1 ^ 2 for x — > 0. We unraveled the following behavior 
(Fig©: 

i) (r 2 ) 1 / 2 diverges for x —> 0,1 in [3 

ii) (r 2 ) 1/2 -> for x -> 1 in jF 
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Hi) (r 2 ) 1 / 2 decreases from the value of ~ 5 GeV 1 (r; 1 fm), at x — > 0, to ~ 2.5 GeV 1 
(« 0.5 fm), at x —* 1 in our model. 

From the behavior of the rms radius with x, one can surmise what components in b space 
the form factor is dominated by. While (r 2 ) 1 / 2 at x > 0.5 is entirely given by b < 2.5 GeV^ 1 
(0.5 fm) for models characterized by a hard component (ii) and Hi)), in case i) saturation 
is not reached (Fig. 10), that is the rms radius is large even at large x. On the other side, 
by studying the average value of x, x ave , calculated using the distribution H vs. A fFigslTTI 
and 112)) . one finds that the hard type functions are not dominated by large x at large A, 
rather all x components seem to contribute (Figll2|). 

In summary, although it is becoming clear that the physics governing the hadronic form 
factors, and more generally exclusive processes at large momentum transfer, is given by a 
non-trivial blend of soft and hard components, current descriptions using GPDs are not 
completely satisfactory. 

The next question is how one would be able to test the transverse spatial structure of 
hadrons, given the fact that GPDs are quite elusive objects from the experimental point of 
view. A potentially powerful means to single out the size of hadronic configurations could 
be given by combining CT/nuclear filtering studies with GPDs as suggested in this paper. 
Based on the so far little explored concepts of IPPDFs and GPDs, we can now test the 
effect of nuclear filtering in an exclusive reaction by building a "filter" in the definition of 
the transparency, Ta- This procedure eliminates some of the model dependence implicit in 
many predictions in a somewhat similar way to what proposed e.g. in Ref . . It however 
also provides a way of explicitly extracting the transverse radial dependence of the hadron- 
hadron cross section in the nuclear medium. While an analysis extended to all available 
data will be carried out in a upcoming paper j24|], we have presented here a few hypothetical 
scenarios obtained with our procedure. 

In addition to the main results of this paper, interesting new developments will concern 
the extension of our calculations to the neutron and the polarized CcLSG, clS well as the 
exploration of the skewedness, £, dependence. Most importantly, the use of nuclei and CT 
experiments might unravel an alternative method to measure GPDs by evaluating nuclear 
dependent contributions beyond the standard assumption of factorization into a nucleon 
times a nuclear part for the nuclear scattering amplitude. Such terms can be shown to 
be directly proportional to GPDs and a study examining their relative importance as A- 
dependent contributions is on its way |24 |. Both the studies initiated in this paper and the 
suggested future studies could be tested ideally both at Jefferson Lab j^| including its 12 
GeV upgrade, and at newly planned Electron Ion Colliders with high luminosity and energy 

0- 

In conclusion, studies like the one presented here using the new concept of GPDs, will 
both improve our knowledge of nuclear filtering phenomena and allow for a more detailed 
understanding of the transverse components involved at large momentum transfer. In par- 
ticular, we hope to have provided a connection between b— and k— space that will help 
to systematically address both the role of Sudakov effects in the endpoints of the hadronic 
wave function [5J, 1 1(1 . and the role of power corrections in the large longitudinal momentum 
regions hinted in [37| . 
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APPENDIX A 

The coefficients defining H(x, A) in Eq. (J26j) are given by: 

pi = 2.780 a; a2 (l -x) 4 - 2 (l + 0.5 x - 3 ) (Ala) 

f 0.5695 - 0.1896 x + 0.6885 x 2 x < 0.4 

p2 = < ~ (Alb) 

F \ 0.605 x > 0.4 v ' 

p3 = 1.118x- a098 (l-x) 2 - 24 (Ale) 

f x < 0.87 , A1J . 

pA = < ~ (Aid) 

F 1 -0.1076x + 0.1079 x > 0.87 v ; 
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FIG. 1: (a) Amplitude for DVCS at leading order in Q 2 ; (b) The same amplitude in a two compo- 
nent model. 
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FIG. 2: The transverse distance between the active quark (upper blob) and the center of momentum 
of the spectator quarks (lower blob) - y; between the active quark and the center of momentum 
of the hadron (full circle) - b, and between the spectator quarks and the center of mass of the 
hadron - hx/(l — x) = y — b. r% and T2 are the space vectors for the active quark and the center of 
momentum of the spectator quarks, respectively, taken in the system of coordinates (x,y,z). For 
simplicity, b is drawn parallel to the x-axis. Assuming IPPDFs that yield values of (b 2 ) decreasing 
as x — > 1, there are three possible behaviors for y and r (Ea. (fTT)|0 : (a) b oc (1 — x) 7 ,7 > 1 
y, r decrease as x increases; (b) b oc (1 — x) 7 , 7 = 1 => y, r are constants independent of x; (c) 
b oc (1 — x) 7 , 7 < 1 => y, r increase as x increases. 
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FIG. 3: The function | <£(k) | 2 , defined in Ea. ()24j) . Full line: model calculation discussed in 
this paper; dotted line: the soft wave function introduced in Ref . [33] : dashed line: model of 
Ref.f.Sfl. including a hard component in $; dot-dashed line, the "semi-hard" distribution proposed 
in RefH- 
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FIG. 4: The GPD, H(x, A), calculated according to Eas. l(22a|) and (J2HJ), plotted vs. A for different 
values of x. Typical values, x = 0.07, 0.36, 0.88, are shown by the arrows. 
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FIG. 5: The form factor, -Fi multiplied by A 4 , calculated in three different models: full line, this 
paper; dashed line, 0; dot-dashed line, 0. Results are compared with the dipole form factor 
(triangles). 



22 



1.2 



^1.1 

a 
i 



< 



0.9 



0.8 



0.7 



0.6 



0.5 



k max = 0.7 GeV 
k max = 1 GeV 
k = 2 GeV 



a gaussian (k =0.7 GeV) 




▲ AAA 







10 



15 



20 



25 



A 2 (GeV 2 ) 



FIG. 6: Dominance of different k components in i*i(A 2 ), showed using the ratio 
-Fi(A 2 , k max )/ Fi(A 2 ). The numerator was obtained by setting the upper limit of integration in 
Eq . (|22c|) to different values of k = k ma x- The dot-dashed line is obtained for k ma x = 0.7 GeV; 
the full line corresponds to k max = 1 GeV, and the dashed line to k m ax = 2 GeV. The triangles, 
obtained by setting k rnax = 0.7 in the model of Ref . [sH] , clearly show the saturation of the ratio at 
low values of k. 
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FIG. 7: IPPDF, q(x,b) obtained by two-dimensional Fourier transforming Eg. 1)26(1 . plotted vs. b. 
The curves shown in the figure correspond to different values of x : 0.01,0.45,0.8,0.9. 
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FIG. 8: Modeling the proton radius, (a) (6 2 ) 1 / 2 in GeV -1 , Eq.©, plotted vs. x, calculated 
in different models: full line, this paper; dotted line 3j|; the dashed and dot-dashed lines were 
calculated using the parametrization proposed in [22 1. This parametrization describes the A 2 



m 



behavior using a gaussian as in [35 h owever the argument which is oc (1 — x)A in [35|| is suggested 

37| . The dashed line and the dot-dashed line were obtained by 



to be (1 - x) n A 2 , with n > 2 

setting n = 2 and n = 3, respectively, (b) Transverse radial quantities: full line, y = (y 2 ) 
Eq.@; dots, b = (b 2 ) 1 / 2 , Eq.©, dashes 
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Eg . ([TO]) . All quantities are calculated in 



the diquark model. For comparison r calculated in the model of |35j|, based on a "soft" type wave 
function, is also shown. 
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the dotted line to b max = 1 GeV -1 ; the dot-dashed line to b max = 0.5 GeV" 1 . Lower panel: ratio 
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b. The continuous lines were obtained using the function in the upper panel. The big dots were 
obtained using the model of Ref . jsjj • Notice that the ratio saturates much faster as a function of 
x for the "hard" distribution adopted in this paper, than for the model of [35I ] . 
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FIG. 11: Average value of x, x ave , calculated using different wave functions in H(x, A) as a funct 
of A. The full line corresponds to the spectator model; the dashed line to the model of Ref.f 



the dotted line to Ref.LJal, the dot-dashed line to Ref. 



37|, for n = 2 (see also Fig®. 
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A (GeV) 

FIG. 12: Comparison of different x components contributions to H (x, A) in the spectator model 
(full lines), and in the model of Ref.|35| (dashed lines). The latter, although being soft, shows the 
dominance of (small) large x components at (small) large A. In the spectator model, regardless 
of its hard k components, also small values of x contribute to large A. The arrows in the figure 
indicate the three different values of x for which the curves were calculated. 
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FIG. 13: Possible scenarios for the effect of nuclear filtering and the onset of CT based on a 
simplified model described in the text. The ratio R of GPDs in a nucleus and in free space is 
plotted vs. Q 2 = A 2 , for two different values of x: x = 0.9 , up per panel and x = 0.3, lower 
panel. For each panel we consider: (a) a oc (1 — x)/x, as in Ref . |35l] (dot-dashed and dashed lines); 



(b) a oc (1 — x) 2 , as in Ref. 
correspond to a large value of the filter's size 
bmax- The effect of the latter is to produce a su] 
at large x for distributions with a large radius |3 



371 ] (full and dotted lines). Moreover, the full and dot-dashed line 
the dashed and dotted lines to a small value of 
ression in the ratio, which is more pronounced 



30 



